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1.  Introduction 


Our  interest  is  in  computing  equilibria  for  the  Walrasian  general 
equilibrium  market  where  there  are  a large  number  of  goods  (say,  two 
hundred)  and  a small  number  of  traders  (say,  five).  The  production  and 
consumption  sets  are  assumed  to  be  closed  polyhedral  convex  sets  and 
utilities  are  assumed  to  be  piecewise  linear  and  concave;  the  algorithm 
is  designed  to  take  as  much  advantage  of  this  "convex  piecewise  linear 
structure"  as  our  knowledge  enables.  More  general  models  of  trade  with 
production  and  consumption  sets  that  are  merely  closed  and  convex  and 
with  convex  preference  orderings  in  lieu  of  concave  utilities  can  be 
approximated  very  closely  by  our  "convex  piecewise  linear  model"  of 
trade,  see,  for  example,  Debreu  [7],  Kannai  [16],  and  Mas-Colell  [20]. 

Our  approach  to  the  problem  is  motivated  from  a number  of  sources, 
namely,  the  attempt  of  Mantel  [19]  to  solve  the  "convex  piecewise  linear" 
model  of  trade,  the  fixed  point  methods  of  Scarf  [24]  for  computing 
economic  equilibria,  a model  of  Mas-Colell's  in  [9]  which  shows  that  the 
equilibria  cannot,  in  general,  be  computed  in  a finite  number  of  steps, 
the  existence  proofs  as  those  of  Arrow  and  Hahn  [l],  Negishi  [2l]  and 
Shapley  [25]  which  are  based  upon  the  household  simplex  rather  than  the 
price  simplex,  the  existence  proof  of  Kreps  [17]  which  uses  the  labeling 
theorem  of  [8]  and  labels  based  upon  payments,  and  last  but  not  least,  the 
computational  results  of  Ginsburgh  and  Waelbroeck  [12,13,14] . The  later 
reported  surprising  convergence  rates  by  using  the  simplex  method  together  with, 
in  effect,  manual  adjustments  of  the  weights  on  utilities  to  obtain  an 
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approximate  equilibrium  via  welfare  optima.  In  a similar  manner  our 


algorithm  uses  the  simplex  method,  but  uses  new  fixed  point  methods, 
namely  [8],  to  adjust  these  weights  automatically  and  to  thereby  compute 


an  approximate  equilibrium  or  show  that  the  traders  are  not  resource 
related. 


The  model  of  trade  we  have  adopted  encompasses  in  a natural  way 


negative  or  zero  prices  and  bankruptcy.  Our  treatment  of  bankruptcy 
is  in  the  sense  of  Arrow  and  Hahn  [1].  Our  development  will  utilize 
the  notion  of  quasi-equilibrium  or  compensated  equilibrium  as  in 
Debreu  [5]  and  Arrow  and  Hahn  [1]. 

The  authors  would  like  to  express  indebtedness  to  Gerard  Debreu 
and  Robert  Wilson  for  helpful  discussions.  Wilson  [26]  has  recently 
announced  an  approach  for  solving  the  "convex  piecewise  linear"  model  of 
trade;  our  approach  is  quite  different  from  his  especially  in  view  of 
the  fact  that  the  methods  were  conceived  over  the  same  period  amidst 
considerable  interaction. 

It  is  assumed  that  the  reader  is  familiar  with  the  simplex  method 
and  the  duality  theory  of  linear  programming,  see  Dantzig  [3]. 


2. 


Traders.  Goods,  and  the  Market 


Consider  a market  with  n traders  and  m goods;  let  i and  j 
in  V = [1,  n]  and  g in  n ^ [ 1,  . . . , m]  index  the  traders  and 

goods,  respectively.  We  note  that  the  use  of  the  word  goods  is  in  a very 
broad  sense,  see,  for  example,  Debreu  [4], 

Each  trader  is  endowed  with  a bundle  of  voods  b^  (a  1 x m vector) 
a sociology,  technology,  and  productive  capability  A^  (an  m x h^  matrix) 
and  a utility  c^  (a  1 x h^  vector).  Without  access  to  the  market, 
each  trader  is  modeled  as  a linear  program 


(1) 


max:  c^x^ 

Xi 

s/t:  A1xi  = bt 


x1  > ° 


This  linear  program  may  be  infeasible  reflecting,  for  example,  that 
domestic  agriculture  may  be  insufficient  to  feed  the  population.  By 
convention  we  shall  take  the  gth  component  of  b.,  b.  , to  be  positive, 

4 ^ O 

zero,  or  negative  according  to  whether  the  trader  has  a positive,  zero, 
or  negative  inventory  of  the  gth  good. 

A price  n for  good  g may  be  positive,  zero,  or  negative. 
Given  access  to  goods  of  the  market  at  prices  it  = (jt^,  ...,  « ) trader 
i is  modeled  as  the  linear  program: 


5 


4—  " a I n iM*l  i I 11  "fc  ■ 


max:  c^x^ 

Xi 

(2)  s/t:  = bt  + yt 

nyi  <0  xt  > 0 

Again  the  trader  attempts  to  maximize  his  utility,  but  now  he  is  permitted 
to  buy  and  sell  in  the  market  to  enhance  his  effort.  The  vector  y^  is 
the  vector  of  net  trades  of  trader  i with  the  rest  of  the  market,  and  the 
budget  constraint  < 0 requires  that  he  not  spend  more  for  goods  than 

he  earns  in  the  market.  The  program  (2)  may  be  infeasible  or  even  yield 
utility  unbounded  above. 

For  a given  n and  x^  we  shall  refer  to  ny^  = ^(A^x.  - b^) 
as  the  (net)  payments  of  trader  i.  According  to  whether  the  payments 
are  nonpositive  (iryi  < 0),  zero  (jty^  = 0),  or  nonnegative  (ny^^  > 0) 
we  shall  say  that  trader  i is  within  his  budget,  has  balanced  his  budget, 
or  is  bankrupt. 

An  equivalent  statement  of  (2)  is 
max:  cixi 

(3)  Xi 

s/t:  itA.x,  < xb.  x.  > 0 

l i — i i — 

Clearly  x^  solves  (3)»  if  and  only  if  (x^,  y^)  solves  (2)  where 
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So  far  it  appears  that  our  model  permits  only  linear  utility ^ but 
we  show  in  Section  6 that  convex  piecewise  linear  utilities  are  a special 
case.  Also  in  Section  6 we  show  that  production  is  included  even  though 
it  does  not  appear  in  our  model  explicitly. 


’ T-  * 


3.  Quasi- Egai libria 

Let  us  define  our  notion  of  a solution  to  the  market,  namely,  a 
quasi-equilibrium  (with  bankruptcy)  ( see, for  example,  the  compensated 
equilibrium  of  Arrow  and  Hahn  [1]).  By  a partition  a u P of  the 
trader  set  v we  mean  that  a (J  P = V,  a Cl  P = p,  a / p,  and  p may 


or  may  not  be  empty. 


x * (*!, 


• . x ) 
* n ' 


We  define  (5 :,  P,  x,  jr)  to  be  a quasi- equilibrium  if 

a)  a U P partitions  the  trader  set, 

b)  x is  feasible  that  is,  £ A.x.  = £ b.  and  x.  > 0. 

c)  each  O.  trader  has  nonpositive  payments  and  subject  to  these 
payments  he  has  maximized  his  utility,  that  is,  itfA^  - b^  < 0,  and 

Xj.  maximizes  cixi  subject  to  i^A^  - t^)  < ^(A^  - b^)  with 
x.  >0,  and 

d)  Each  P trader  has  nonnegative  payments  and  has  minimized 

his  payments,  that  is,  i^Aj^  - bi)  > 0 and  5L  minimizes  it^x^  - b^ 

subject  to  x >0. 

i — 

Observe  that  both  (a,  p,  x,  it)  and  (a',  p',  x,  it)  could  be 
quasi-equilibria  with  0!  ^ Q;',"  the  point  being  that  minimizing  payments 
may  have  by  coincidence  maximized  the  utility  subject  to  the  budget  or 
vice  versa. 

If  (a,  p,  x,  jt)  is  a quasi-equilibrium,  but  (v,  p,  x,  it)  is 
not,  then  we  call  (0!,  p,  x,  n)  a proper  quasi-equilibrium.  On  the 
other  hand  if  (v,  p,  x,  jt)  is  a quasi-equilibrium  then  we  call  ( x,  jt) 
an  equilibrium;  for  emphasis,  note  that  (x,  it)  is  an  equilibrium,  if 
and  only  if  x is  feasible  and  each  trader  with  >£  has  maximized  his 
utility  subject  to  his  budget. 
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Fundamental  questions  are: 


i)  When  does  an  equilibrium  exist? 

ii)  How  does  one  compute  an  equilibrium  or  show  that  it  does  not  exist? 

Current  theory  does  not  supply  complete  answers,  but  we  shall  cover  the 
major  developments  known  to  us. 

In  general  not  even  a quasi-equilibrium  exists  much  less  an  equilibrium. 
The  next  three  conditions  guarantee  the  existence  of  a quasi-equilibrium. 
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I)  Feasibility:  We  say  the  market  is  feasible,  if  there  is  an 

x > 0 such  that  E A.x.  = E b . . □ 

— V i i Vi 


II)  Nonsatiation:  We  say  that  the  market  possess  the  nonsatiation 

property,  if  each  c^  contains  at  least  one  positive  element.  □ 


III)  Finite  Utility:  We  say  that  the  market  has  finite  utility, 

if  for  each  i the  program 

max : c . x . 

1 l 

s/t:  E A x = E b x > 0 

V J \i  -* 

has  a finite  objective  value.  □ 

If  a quasi-equilibrium  is  to  exist,  clearly  the  market  must  be 
feasible;  the  other  two  conditions,  though  not  necessary,  also  seem 
inoccuous.  Each  of  these  three  conditions  are  easily  checked. 


Assumption:  Throughout  the  remainder  of  the  paper  we  assume  that 


the  market  is  feasibile  (I),  possess  the  nonsatiation  property  (II),  and 


has  finite  utility  (III).  □ 


In  Section  5,  we  supply  a constructive  proof,  in  the  limiting  sense, 
to  the  following  theorem. 


Theorem  1:  A quasi- equilibrium  exists.  □ 


Observe  that  the  price  it  of  a quasi-equilibrium  is  nonzero;  this 


follows  from  the  nonsatiation  assumption  and  the  fact  that  OC  is  nonempty. 
Suppose  that  trader  i is  individually  feasible,  that  is,  ( 1)  is  feasible, 
and  that  (a,p,x,jr)  is  a quasi-equilibrium;  then  trader  i is  not  bankrupt, 


that  is,  trader  i is  within  his  budget.  Otherwise  nb^  < 0 and  hence, 
with  jtA^  >0  we  see  that  ( 1)  is  not  feasible.  In  general  one  cannot 


conclude  that  an  individually  feasible  trader  is  of  CL  type. 


The  next  lemma  shows,  in  particular,  that  a quasi-equilibrium 


without  bankruptcy  is  O-Pareto  efficient  and  that  an  equilibrium  is  in 


the  core.  In  a quasi-equilibrium  observe  that  an  Ct-trader  cannot  achieve 


the  same  utility  at  a lower  payment  due  to  the  nonsatiation  assumption  (II); 


also  observe  that  trader  i can  be  a 0- trader,  if  and  only  if  ,'A^  > 0, 

irA.x.  = 0.  and  rtb.  < 0. 
i l » l — 


Lemma  2:  If  (a,  p,  x,  n)  is  a quasi-equilibrium,  then  there 


does  not  exist  an  CL  c OLf  p c p^  x^  > 0 for  i in  OL  U p,  and  6.  > 1 


for  i in  p such  that 


■■■■■ 
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* Z (A.x  - b ) > 0 
oue 


z 

a 


A.x. 
1 1 


f£5.Ax.  = y b . + E S.b. 
p 1 1 1 a 1 p 1 1 


c.x.  < c.x.  for  all  i in  a 

l l — ii 

and 


c.x.  < c.x.  for  some  i in  0! 

ii  ii 


Proof.  If  so.  then  n(A.x.  - b.)  >jt(A.x.  - b.)  for  all  i in 
* ii  i'  — v i l i' 

0!  and  with  a strict  relation  for  some  i in  OC.  Since  6.  Jr( A x - b 1 

l ' i i i ' 

> jt(A.x.  - b.)  for  each  6-trader  we  have 

— ' i i i' 


*(£  (A  x - b ) + Z MAx  - b ))  > 0 

a p 1 


which  contradicts  the  feasibility  of  for  i in  OL  U |3.  □ 


The  existence  of  a proper  quasi-equilibrum  might  lead  one  to 
believe  that  an  equilibrium  does  not  exist.  For  consider,  the  Cfc- traders 
are  within  their  budgets,  and  the  (3- traders  have  foresaken  their  utilities 
and  focused  on  balancing  their  budgets  which  they  are  just  able  to  do, 
if  at  all.  The  following  example  illustrates  that  the  model  may  have 
both  a proper  quasi-equilibrium  and  an  equilibrium. 
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cL  = (1,  0)  c2  = (1) 


One  can  check  that  x^  = ( 1^  1)  f x,^  = 1 and  n = ( lf  -1)  is  an  equilibrium. 
However  a = { 1} } £ = (2),  = (2}  0)  f Xg  = 0,  and  it  = (1,  0)  is  a 

proper  quasi-equilibrium^  the  (3- trader  has  minimized  his  payments  and 
has  thereby  managed  to  balance  his  budget;  however^  he  has  not  maximized 
his  utility  subject  to  his  budget. 

Computation  of  an  approximate  quasi-equilibrium  in  a finite  number 
of  steps  involves  some  technical  difficulties  that  we  have  been  unable  to 
surmount.  We  are  thus  lead  to  the  "indirectly  resource  related"  notion 
of  Arrow  and  Hahn  [l], 

(IV)  Resource  Related:  The  traders  of  the  market  are  defined  to 

be  (indirectly)  resource  related  if  no  subset  of  traders  is  satiated  by 
the  goods  owned  and  produced  by  the  remaining  traders.  More  precisely ^ 
we  say  the  traders  are  resource  related  if  for  any  x with 


and  any  partitions 
&^  > 1 for  i £ p 


a U P of  the  traders  with  p / 0,  there  exists  numbers 
and  x > 0 such  that 


i)  Z Ax  + Z 6,A.x.  = Z b.  + Z &.b.  , 
a 11  p 111  a 1 p 11 


ii)  cixi  < cixi  for  311  i in  a. 


iii)  cixi  < cixi  for  some  i in  a.  □ 


Theorem  3:  If  the  traders  are  resource  related,  then  an 

equilibrium  exists. 

Proof : There  is  a quasi-equilibrium  [0Lf  B,  xf  it)  by 

Theorem  1.  If  p / 0,  then  according  to  Lemma  2,  the  resource  related 
condition  fails  for  (af  B,  x) . □ 

The  algorithm  described  in  Section  7 will  compute  (in  a finite 
number  of  steps)  an  approximate  equilibrium  or  show  that  the  resource 
related  condition  (IV)  fails. 
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4.  The  Aggregate  Program  with  Weighted  Objective 


For  our  proof  and  algorithms  we  consider^  in  the  spirit  of  welfare 
economics^  the  aggregated  production  capability  of  the  market  together 
with  an  objective  function  obtained  by  weighting  the  traders 
utilities. 

Let  S be  the  ( n- 1) -simplex  {(0^,  0R)  >0  : Zy  6^  = 1). 

For  0 £ S consider  the  linear  program: 

f(0)  = max:  Z 0.c.x. 

' ' ill 

x v 

(4)  s/t:  Z A x£  = Z b x > 0 

V v 


We  observe  that  f is  a finite  continuous  piecewise  linear  convex  function 
on  S. 

The  dual  of  (4)  is 

min:  jt(Z  b ) 

jt  v 

(5) 

s/t:  »t  A , > 0 . c . i£v 

i — i l 

For  any  0 £ S we  know  that  the  optimal  objective  values  of  (4)  and  (5) 
are  equal  and  finite. 
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We  define  three  point  to  set  maps,  namely,  X : S -+R^-hi  , jj  ; s ->  Rm 

and  P : S ->  R . For  each  6 £ Sf  X(6)  is  defined  to  be  the  set  of  optimal 

solutions  to  (i+),  Jl(e)  is  defined  to  be  the  set  of  optimal  solutions  to 

(5)>  and  P(0)  is  defined  to  be  the  set  of  payment  vectors  p(x,it,0) 


p(x,rr,e)  = 


jr(AlXl  - bL) 


jt(A  x - b ) 
n n n' 


where  x £ X(  9)  and  n € Il(  0) . We  observe  that  6 c x = nA  x for 

i i i i i 

any  i £ V,  x £ X(  Q)  f and  jt  £ n(  0)  . Hence  p(X>n,0)  is  also  of  form 


01C1X1  " nbl 


Bex  - jtb 
n n n n 


for  x € X(6)  and  * £ n(0).  If  P = ( Pp  pn)  £ P(6),  then  clearly 


Z p.  = o 

V 


Lemma  4:  Assume  x and  it  are  in  X(0)  and  11(6).  If 

P1(x,tt,6)  £0  for  0i  > 0 and  p^x.ir.O)  _>  0 for  0 = 0,  then 
(a,  B,  x,  it)  is  a quasi-equilibrium  where  a = (i  : 0.  > 0}  and 

B>  = U : e1  = o). 
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Proof:  Obviously  tt  U P is  a partition  of  v and  x is  feasible. 

For  i in  5 we  have  ivA.  > 0. c.  and  uA.x.  - 0.c.x..  If 

i — ii  ii  ill 

itA.x.  - nb.  < nA.x.  - itb..  then  0,c.x.  = nA.x.  > 0.c.x.  and  hence 

ii  l—  ii  l'  ill  li—iii 

trader  i has  optimized  subject  to  his  payment  level.  For  i in  £ 
we  have  nA.  > 0.  itA.x.  = 0,  and  nb.  <0.  □ 

In  particular,  if  we  can  prove  that  there  is  a Q > 0 for  which 
P( 0)  contains  a zero,  then  we  have  obtained  existence  of  an  equilibrium. 
Note  that  if  q £ P(0),  then  the  task  of  computing  an  x £ X( 0)  and 
it  £ Il(  0)  with  p(x,it,0)  = q is  merely  a matter  of  solving  a linear 
program. 

The  next  theorem  describes  the  continuity  properties  of  X,  H,  and 

k 

P.  A point  to  set  map,  say  X,  is  called  upper- semi- continuous  if  0 

k 

tends  to  0,  x in  X(  0 ) , and  x tends  to  x,  imply  x is  in  X(  0) . 

A point  to  set  map,  say  H,  is  called  lower  semi-continuous,  if  rr  in 

k k k 

H(0)  and  0 tends  to  0 imply  there  exist  it  in  H(  0 ) with 

k 

jt  tending  to  jt.  A map  which  is  both  upper  and  lower  semi-continuous 

is  called  continuous.  We  shall  need  the  following  lemma  which  can  be 

proved  using  the  fact  that  the  objective  of  a linear  program  is  a continuous 

k k 

function  of  its  right  hand  side  over  the  feasible  range;  let  1)  Bz  < d 
where  k = 1,  2,  . . . be  a sequence  of  inequality  systems  in  the  variable 
z. 


Lemma  5:  If  D is  a diagonal  matrix,  Dk  tends  to  D^,  D°°  has 
a positive  diagonal,  dk  tends  to  d°°,  and  DkBz  < dk  has  a solution  for 
k = 2,  . .. , then  D Bz  < d has  a solution,  and  for  any  such  solution 

z there  is  a sequence  zk  tending  to  z'  with  DkBzk  < dk.  □ 


lit 


Lemma  6: 


For  each 


empty,  convex, 
is  continuous. 


and  closed. 


0 the  sets  X(  0)  , Jl(  0) , and  P(  0}  are  non- 
X and  P are  upper  semi-continuous,  and  u 


Proof:  We  shall  only  prove  that  P is  upper  semi-continuous; 

k k 

the  remaining  facts  are  much  easier  to  establish.  Assuming  p G P( 0 ) , 

pk  -»  p°°,  and  0 k ->  0°’  we  need  to  show  that  p°°  £ P (©“).  Select 

xk  € X(0k)  and  nk  € Il( 0k)  so  that  p(xk,  irk,  0k)  = pk.  If  xk  -*  x°° 

and  itk  ->  n°°,  then  x“  £ X(  0 ) , jt°°  t ll(  0 ) and  p(xa’,  jt°°,  G°°)  = p°’,  and 

the  result  is  established.  The  difficulty  arises  because  xk  and  or 
k 

jt  may  not  converge. 

Ic  k 

For  k = 1,  2,  . ..,  oo  the  condition  that  p is  in  P(0  ) is 
equivalent  to  the  requirement  that  (6k)  have  a solution 


(6k) 


(ai) 

k k 

0.c.x.  - itb . = p.  , 
ill  i l * 

i € v 

(b) 

Z AiX.  = Z b 

V V 

x > 0 

(c) 

* 

CD 

II 

JQ 

(d) 

nA.  > 0kc.  , 
i — i l > 

i € v 

Observe  that  by  summing  each  (ai)  and  (c)  we  get  Yy  0^^  = f( & ) 
since  Zy  pk  = 0 for  k = 1,  2,  . ..,  oo. 

Now  observe  that  by  the  assumption  of  finite  utility  we  cannot  have 

k k k k 

(?iCiXi  ->+oo,  and  consequently,  we  cannot  have  O^c^x.  -0°. 


Pick  a subsequence  of  k so  that  each  0*  and  converges. 


k k 

Now  apply  Lemma  5 to  (6k)  where  D = I,  x is  set  equal  to  x , and  where 


only  n is  regarded  as  a variable. 


We  conclude  that  there  is  a sequence  n tending  to  n where 


(x^,  jt^)  solves  (6k)  for  k = 1,  2,  . . . . If  0^  ->  0,  then  clearly 
k k 

^iCiXi  ^ow  apP^  Lemma  5 to  the  system  formed  by  the  (ai^)  with 

0°°  >0  and  (b)  where  n = x is  the  only  variable  and 

J 


We  get  a sequence  x ->  x solving  the  reduced  system  for 

Me  ~k 

k = 1,  2,  . . . , oo.  We  now  have  a solution  (x  , jt  ) to  (6k)  tending  to 
(2s,  rt)  and  the  result  is  complete.  □ 


Existential  and  computational  difficulties  are  caused  by  the  fact 
that  the  X(0),  ll(0),  and  P(  0)  may  be  unbounded  and  we  move  to  define 

/S  /\  /N 

selections  X(0),  11(0),  and  P(0)  from  them.  Given  two  ^-vectors 
and  v we  say  that  u is  lexico  greater  than  v,  if  there  is  an 
i = 1.  ....  £ such  that 


u 


m — m 


and 


U = V 

j J 


ui > vi 


for  1 < j < i_1 


Given  a set  of  u's  the  lexico  maximum  (resp.  minimum)  is  the  u which 
is  greatest  (resp.  least)  in  this  ordering. 

A 

Let  X(0)  be  the  unique  solution  to  the  (lexico)  linear  program 

lexico  max:  (E  “ x) 

s/t : E A^  = Z b,  ( x > 0 . 


Here  one  first  maximizes  E Q^c^x^,  then  minimizes  x^,  then  x^, 

etc.,  and  finally  xnh  . Obviously  X(0)  is  an  element  of  X(  6) . 

n 

Furthermore,  it  is  easily  verified  that  X( d)  uses  linearly  independent 
columns,  that  is,  the  columns  of  (A  ...,  A ) corresponding  to 
positive  elements  of  X( 6)  are  linearly  independent.  It  follows  that 
X(S)  = (X(0)  : 6 e S}  is  finite. 

Select  A an  m X h matrix  with  hrt  > 0 so  that  (A„.  ...  A ) 
u u u “ ' 0*  * n7 

has  rank  m and  so  that 

Vo  + l Vi = l bi 


x0  > o,  x > 0 


implies  Xq  - 0.  Such  an  AQ  is  easily  computed  with  elementary  row 
operations. 
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e™)  and  let  B be 


rr ? 


Let  [ e ] be  the  column  vector  (e,  e* , ...f  e”1)  and  let  B be 
any  m x m submatrix  formed  using  linearly  independent  columns  of 
<V  A n)  • each  Q define  il(0)  to  be  the  unique  solution  to 


* 

the 

dual 

of 

1 

max: 

Z e icixi 

1 

V 

1*. 

s/t: 

Aoxo  + z 

AiXi  = ^ b 

[ > 

V 

V 

x.  > 0, 

x > 0 

i 

0 — > 

1 

IL 

for 

all 

small  e > 0 

where  the  variables  are 

I 

f 

1 

n(e) 

is 

defined  to  be 

the  unique 

solution  to 

program 


lexico  min:  (it(E  b ) } itB) 
V 


s/t:  jtAq  > 0 

«Ai  > 9lci  i € v . 

It  follows  that  n(6)  is  in  n( 6)  and  that  P(0)  is  in  P(e)  where 
we  define  P( Q)  to  be  p(X(e),  R(6) , &) . 

/n  in 

Lemma  7 : The  function  n : S -> R is  continuous  and  nonzero. 
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J 

1 


*r 


Proof.  ll(  6)  £ b.  = £(&)  is  a continuous  function  of  0.  by 

v i ’ 

A 

the  same  reasoning,  n(0)B  is  a continuous  function  of  0,  and  it  follows 
that  n is  continuous.  Since  II ( 6 ) f°r  » H(0)  ^ 0.  □ 


Since  X(S)  and  IL(S)  are  bounded,  P(S)  is  bounded.  In  general 

/-V  /N 

X and  P are  not  continuous. 

Finally,  we  observe  that  X(  0)  and  n(0)  can  be  computed  by 
appropriate  application  of  the  simplex  method  to  the  ( lexico)  linear 
program 


5.  Existences  of  Quasi-Equilibria 

/N 

Using  the  function  P,  the  proper  labeling  theorem  of  [8],  and  Lemma  4 
we  prove  the  existence  of  a quasi-equilibrium. 

Let  l : S -»  Rn  be  a function^  continuous  or  not,  on  the 
(n- 1) -simplex  S.  If  the  origin  of  Rn"^is  in  the  convex  hull  of  1(C)  for  a 
subset  C of  S we  say  that  C is  g-complete.  If  every  neighborhood 
of  0 is  an  ^-complete  set  we  call  8 an  ^-complete  point.  The  proof  for 
the  following  result  is  found  in  [6]  or  [8]j  this  result  can  be  regarded  as 
equivalent  to  the  Brouwer  fixed  point  theorems. 

Lemma  8:  If  the  vertex  set  of  S is  ^-complete  and  no 

(n-2)-face  of  S is  ,0-complete,  then  there  is  an  ^-complete  point.  O 

Our  next  step  is  to  define  a labeling  Z : S -*  { z € Rn  : Y^  z ■ = 0} 
related  to  our  market.  For  8 on  the  boundary  of  S define  Z(  6)  by 


V0)  = 

_1 

CT 


if  0 > 0 


if  ejL  = o 


where  a is  the  number  of  i's  for  which  8 = 0.  For  8 in  the  interior 
of  S define  1(8)  to  be  P(  0) . 

We  see  that  the  vertex  set  of  S is  i-complete  and  that  no 
(n-2)-face  of  S is  ^-complete,  and  we  may  conclude  that  there  is  an 
i-complete  point. 
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Lemma  9 : 


For  each  ^-complete  point  0 there  is  a p t p(g) 
with  Pt  < 0 if  et  > 0 and  Pt  > 0 if  § = 0. 

Proof : Froin  the  definition  of  an  ^-complete  point  it  follows 


that 

there 

ik 

is  a subsequence  of  0 in  S 

ik 

and  A >0 

for  i in  v 

such 

that 

.ik  - , 

0 -»  0 for  each  i in  v as 

k=  1,  2,  ... 

tends  to  « 

and 

Z xlk  «eIk)  - o 

V 


V 

for  k - 1,  2,  ...  . 

By  rearranging  terms  and  selecting  subsequences  if  necessary,  we 


may  assume 

that  0ik  >0  for  i 

in  £ 

, -ik  i 
and  0 f 

0 for 

i in 

where 

C U 

q partitions  v. 

Further  we 

may  assume 

that  A 

ik  • 

->  A 

for  i 

in 

^ik.  i 

V,  f(0  ) -»  p 

for  i 

in 

S,  and  £(Q 

ik.  i 

) -*  q 

for 

q as 

k = 

1,  2,  ...  tends 

to  oo 

where 

and 


We  have 


P1  € p(e) 


< o if  e = 0 , 

i - _ 

9j  > 0 if  > O . 
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Divide  by  s = A and  use  the  convexity  of  P(0) 

witrh 

to  get  p in  P(e) 

Pj  > o 

if  0j  = 0 , 

Pj  < 0 

if  > 0 . 

□ 

Putting  together  Lemmas  4,  8,  and  9 we  have  a proof  of  Theorem  1 
namely,  a quasi-equilibrium  exists.  Lemma  4 also  yields: 

Lemma  10:  If  there 

is  an  equilibrium,  □ 

is  an  i-complete  9 with 

9 > 0f  then  there 

Mm 
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6.  Piecewise  Linear  Utility  and  Production 


Assuming  production  and  consumption  sets  are  polyhedral  convex 
sets  and  that  preferences  are  reflected  in  a concave  piecewise  linear 
utility  function  we  shall  show  that  the  traditional  model  of  competitive 
equilibrium  as  described,  for  example,  by  Debreu  [4]^  is  a special  case 
of  our  model  as  stated  in  Sections  2 and  J.  We  shall  incorporate  pro- 
duction using  a construction  attributed  to  Rader  [22], 

Let  us  suppose  that  there  are  m^  goods  in  the  ( traditional) 

mQ 

market.  Let  the  production  set  be  a set  in  R of  form 

<Vo  : Vo  * V *0  £ °)  • 

Note  that  any  closed  polyhedral  convex  set  can  be  represented  in  this 
manner.  It  is  assumed  that  disposal  is  free 

< 0 implies  q^  is  in  Qq 

that  production  in  the  large  is  not  reversible 

{qg  € Qq  : qQ  > 0)  is  bounded. 

Given  a set  of  prices  itQ  for  the  goods,  the  profit  j^qQ  is  to  be 
maximized  over  qg  in  Qq. 
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Let  us  suppose  that  there  are  n consumers.  We  assume  that  the 

mO 

consumption  set  Q of  consumer  i is  a set  in  It  / of  form 


{E.  z.  : D. z.  = dJ . z >0) 
11  ii  i*  i — 1 


Each  consumer  is  assumed  to  have  an  endowment  of  goods  w.  and  a concave 
piecewise  linear  utility  function  over  the  consumption  set  defined 


U. ( q)  = min  c . . q 
i . _ ii 


J € H, 


where  is  a finite  set.  Observe  that  any  concave  piecewise  linear 

utility  function  can  be  represented  in  this  form.  Further  we  suppose  that 
each  consumer  owns  a fraction  y^  > 0 of  the  production  facility  where 

ri  = u 

We  assume  that  Q.  has  a lower  bound 


q.  in  Q.  implies  q.  > r 
i i l — 


that  there  is  no  satiation  consumption 


^i  ^i  there  is  a 


q.  in  Qi  with  U(qi)  > U( q.) 


and  that  there  is  a q.  in  Q.  such  that  q < w . 

1 l l i 


CJj 


f 


Given  this  s tatement  of  elements  of  the  model  we  now  recast  it 


> 


l 


¥ 1 


* 


into  the  form  of  Sections  2 and  3. 

Given  prices  jt  consumer  i's  problem  is  to 

max: 

s/t:  q.  € Q. 

rtqi  ^ *owi  + ri 

where  it q^  is  consumer  i's  income  from  the  production  and  where 
qg  maximizes  ir^q^  subject  to  q^  € Qq.  This  problem  of  consumer  i 
can  be  restated  as 

max-'  v^ 


V,  < c 

i — : 

Ljqi  ' 

j € 

E.z.  = 

q. 

1 1 

0 

N 

II 

d.  . 

z.  > 0 

i i 

i - 

Ez  = 
0 0 

% 

D z„  = 

d„ 

z > 0 

0 0 

0 ’ 

0 - 

n0qi^ 

rt0wi  + ri’TOqO 

where  the  variables  are  v^,  q^  z z^f  q and  q^.  Observe  that  if 
the  program  has  a solution  (in  view  of  the  nonsatiation  assumption)  or 
if  the  consumer  minimizes  his  payments  then  the  production  facility  is 
operated  optimally. 
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■Mesa*? 


LU  - 


* *■  « *«►  i 1 

•jku i £Wr  * rV 


i 1 


k i 


vi  > °,  vi  > °,  sj  > ° 


Eizi  - riEozo  - wi  = yi 


D z = d 

i i i > 


z.  > 0 


D z = d 
0 0 0 > 


'oi° 


Vi  < 0 


where  the  variables  are  v+,  v^,  s^,  z^,  zQ,  and  y^.  Aggregating 


all  these  variables  except  to  form  the  variable  and  defining 


12  12 

A^,  A.,  b^,  b^,  and  appropriately,  we  can,  again,  restate  the 


consumers  problem  as 


max:  c.x. 

l l 


/ * 1 ,1 

s/t:  A.Xi  = bi  + 


2 2 

AIXi  “ bi 


> 0 


Vi 
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CfFVX'-i,  - ' \L  .*  - ham  - ir 


Now  assuming 


A,  has  m,  rows  we  define  A.  and  b,  to  be  the 
i 1 11 

matrix  and  vector 


< — Rows  1 to  mQ  — > 


0 

A 

0 


2 

i 


0 


i- 1 i 

<* — Rows  ( Yj  m.)fl  to  £ m.  — > 

0 1 0 L 


b2 


l 

0 


Hence,  in  the  model  of  Sections  2 and  5 there  are  m = m. 

’ ~0  l 

goods,  but  only  mQ  of  them  are  traded  at  a quasi-equilibrium  in  view 
of  the  feasibility  requirement. 

Since  there  is  q^  < w^  and  disposal  is  free,  the  market  is 
feasible  (I).  Clearly  c^  has  a positive  component  and  insatiation 
(II)  follows.  Since  each  consumption  set  is  bounded  below  and  the  pro- 
duction is  irreversible  we  obtain  finite  utility  (III).  Thus  the  model 
has  a quasi-equilibrium. 

At  a quasi-equilibrium  obviously  the  prices  jt0  are  nonnegative 
and  nonzero.  Since  each  trader  has  a q^  < w^  the  trader  can  easily  im- 
plement a surplus  budget.  Thus  we  have,  in  fact,  an  equilibrium.  By  making 
assumptions  as,  for  example,  in  Bergstrom  [2]  or  Hart  and  Kuhn  [15]^ 
equilibria  would  not  necessarily  yield  nonnegative  prices. 


7 . The  Algorithm  and  e-Equilibria 


Assuming  feasibility  (I),  insatiation  (II),  and  finite  utility  (III) 
we  have  shown  existence  of  a quasi-equilibrium,  and  our  attention  is  now 
turned  to  computing.  In  general  one  cannot  compute  in  a finite  number 
of  steps  (additions,  multiplications,  and  branches  on  signs)  a quasi- 
equilibrium, for  there  are  markets  with  rational  data  (A,  b,  c) 
where  all  quasi-equilibria  have  irrational  components,  see  [9]  and  Section  8. 
However,  if  there  are  only  two  traders  (n  = 2),  then  a quasi-equilibrium 
can  be  computed  in  a finite  number  of  steps,  see  Gale  [10].  The 

next  logical  attempt  would  be  to  compute  in  a finite  number  of  steps  an 
approximate  quasi-equilibrium.  However  we  have  not  been  able  to  devise 
an  algorithm  and  a satisfactory  notion  of  an  approximate  quasi-equilibrium 
where  any  degree  of  efficiency  could  be  anticipated.  We  proceed  to 
describe  an  approximate  equilibrium  and  an  algorithm  which,  in  a finite 
number  of  steps,  computes  an  approximate  equilibrium  or  computes 
(<2,  P,  x)  showing  that  the  traders  are  not  resource  related. 

Our  notion  of  an  approximate  equilibrium  is  that  of  an  e-equilibrium. 
Towards  normalizing  the  prices  we  let  | it  | = I™  |it  | • Observe  that  if 
(x,  it)  is  an  e-equilibrium,  then  so  is  (x,  gott)  for  any  w > 0 and,  in 
particular,  for  w = |ir|  1 . We  call  (x,  tr)  an  e-equilibrium  if 

a)  x is  feasible,  that  is  A^x^  = Zy 

b)  Each  trader  is  within  e of  balancing  his  budget,  that  is, 
nA^x^  < rtb^  + 1 7T  | e . 

c)  Each  trader  is  within  e of  optimizing  subject  to  his  budget, 

that  is.  c.x.  + e exceeds  c.x.  for  all  x.  > 0 with  nA.x.  < nb.. 

' i i it  i—  i i — i 
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Our  algorithm  consists  of  three  parts,  a scheme  for  generating 
^-complete  sets,  a termination  rule  which  detects  an  e-equilibrium, 
and  a termination  rule  which  detects  that  the  traders  are  not  resource 
related.  The  scheme  for  generating  the  ^.-complete  sets  is  prescribed  in 

[6]  or  [8]  and  we  shall  only  sketch  it  here.  Using  the  labeling  of 

Section  5,  the  scheme  generates  a sequence  of  ^-complete  sets  0*,  ©^  ... 

where  ©k  = (0lk,  ...,  0nkj  c S,  the  diameter  of  the  ©k,s  tend  to  zero 

k+1  k 

as  k tends  to  infinity,  and  © ~ © contains  one  element  for  all  k. 

Given  ©k  at  iteration  k a 0k  is  computed,  a label  £(0k) 

k+ 1 

is  computed,  a pivot  is  made  as  in  the  simplex  method,  and  © is 

thus  determined.  Thus  if  the  9 is  interior  to  S,  the  aggregate 

linear  program  (4)  must  be  solved  to  compute  Jt(0  ).  The  presumption  is 
k 

that  after  the  © 's  get  small  the  solution  of  the  aggregate  program 

for  9 is  easily  obtained  by  using  the  simplex  method  to  post  optimize 
k- 1 

from  9 . Also  we  note  that  the  aggregate  linear  program  will  most 

likely  have  much  special  structure  thereby  lending  it  to  the  vast  store 
of  large  scale  linear  programming  techniques,  see  Geoffrion  [11]  and 
Lasdon  [18]. 

If  ©k  is  interior  to  S,  then  letting  x^k  equal  X(6^k)  and 
ir^k  equal  n(6^k)  we  have 


(8) 


£ UjkA.  xf.xjkb.]  v>k=0, 
j€v  11  1 


where  A^k  > 0 and  £ AJR  = 1.  Define  jt"’  = £ jtjr‘  AJ^,  x"-  = £ x 


Jk 


Jk  -Jk  k 


Jk  %Jk 


and  ek  = ^ ejk  Ajk. 


We  shall  introduce  two  termination  rules.  The  first  detects  an 
e-equilibrium  and  the  second  detects  that  the  traders  are  not  resource 

re lated. 
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Termination  Rule  1:  For  a given  e > 0 if  © is  interior  to 


S and  both 


(*Jk  - 


eJ)c 


Xjk 

i i 


e: 


and 


, Ik  k.  jk 
(ttj  - it  ) A±  xJ± 


have  an  absolute  value  of  e or  less  for  all  i and  j in  v,  then 
terminate  the  algorithm.  □ 


Lemma  11:  If  the  algorithm  terminates  with  Rule  1 at  step  k, 

“ • k k 

(x,  it)  = (x  f n ) is  an  e-equilibrium. 

Proof:  Since  Zv  x|k  = Zv  b.  we  have  A.i^  = 7 b.  and 

we  see  that  x is  feasible.  Next  it  must  be  shown  that  each  trader  is 
within  e of  meeting  his  budget  and  within  e of  optimizing. 

From  ( 8)  we  have 

Z [(it  + (n^k  - it))  A x-?k  - n^k  b ] ?^k  = 0 

j i 1 i 

or 

ifAi*i  = »rb1  + Z («jk  - it)  x|k  Ajk 

or  finally 

itA^x^  < nb^  + | tt  | e 

that  is,  each  trader  is  within  e of  meeting  his  budget. 


then 


From  (8)  and  t^k  c^  x|k  = jt^k  Ai  x|k  and  letting  0 = C k we  have 


or 


E [(©,_  + (e|k  - 5t)  c±  xJtk  - *jk  bt]  *jk  = 0 


©i  = ^ + E (©|k  - c.x|k  Xjk 


or  finally 


*bi 

c x > e . 

'6i 

Ik  i k . . 

We  have  n c^  and  hence  jrA^  > Thus  at  prices  jf  for 

a trader  to  meet  his  budget  we  have  JtA^x^  < itb^,  hence  £Jicixi  < nb^  or 


”bi 

cixi  — j • 

h 


Therefore,  given  his  budget,  trader  i is  within  e of  optimizing.  □ 

■ lr 

We  regard  6 as  a cluster  point  of  the  © if  every  neighborhood 
of  G contains  infinitely  many  of  the  @k.  Of  course,  the  sequence  ©k 
has  at  least  one  cluster  point. 

Lemma  12:  Let  e be  given  and  assume  that  the  sequence  of 

k 

^-complete  sets  © has  at  least  one  cluster  point  interior  to  S. 

Then  the  algorithm  will  terminate  after  a finite  number  of  steps. 
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J - * 


IT 


Proof: 


A 


k occurs  for  which  Rule  1 will  terminate  the  algorithm, 
(if  Rule  2 has  not  already)  since  the  x^'s  come  from  a finite  set  X(S), 
H is  continuous  and  nonzero,  the  diameter  of  6^  tends  to  zero,  and 

t, 

© has  a cluster  point  interior  to  S.  □ 

For  the  second  termination  rule  we  say  that  ©^  is  within  e of 
the  p-face  if  < e for  all  i in  3 and  j in  v.  Let  K be  an 

infinite  subsequence.  For  0 in  S let  us  define  y(  0)  to  be  the  set 
of  (i,j)  such  that 

n(0)  Au  = 0t  ct.  . 

Termination  Rule  2:  Let  & and  the  subsequence  K be  given. 

If  k is  in  K,  solve  the  following  linear  program  for  all  (3,  y(0)) 

k - k 

such  that  © is  within  & of  the  6- face  and  0 is  in  © . 


max:  z 


s/t: 

”Au  * 9icti 

(M)  © r(e) 

"Aij  - eiclJ 

(i,J)  t r(3) 

«i-° 

i € 3 

nbi  < 0 

t € 3 

e±  > z 

i t & 
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The  variables  are  zf  Qf  and  jt.  If  the  optimal  z is  positive^  terminate 
the  algorithm.  □ 

Lemma  13:  If  the  algorithm  terminates  with  Rule  the  traders 

are  not  resource  related. 


Proof : If  the  system  has  a solution  then  we  have  a 0 in  the 
boundaiy  of  S and  n in  n( 0)  where  0 = 0 implies  nb^  < 0.  Hence 
the  following  lemma  applies.  □ 

Lemma  14:  For  fixed  0 in  the  boundary  of  S and  jt  in  li(  0) 
suppose  that  0^  = 0 implies  nb^  < 0.  Then  the  traders  are  not  resource 
related. 


Proof:  Let  a = { i : ^ > 0) , 3 = {i  : 0 = 0},  and  x € X(0). 

If  > 1 for  i in  p and  x > 0 satisfies 


Z Ax  + z 8tA  x -Zb  tZ 


a 


we  have 
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E 

a 


0 . c .x.  < Y jtA.x  . 

1 1 1 " a 11 


= E Jtb  + E &,Jtb  - E & .nA  x 

a e 3 

< E rtb  + E trb. 

a 3 

= E nA  X = E 0.C.X  . D 

a a 

Theorem  15;  Let  e > 0,  5 > 0,  and  K be  given.  Then  the 
algorithms  terminates  after  a finite  number  of  steps  with  an  e-equilibrium 
or  a (a,  3,  x)  showing  the  traders  are  not  resource  related. 

Proof:  In  view  of  Rule  1 we  can  assume  that  0^  does  not  have 

any  cluster  points  interior  to  S.  Clearly  y(S)  is  a finite  set  and 

©(q)  = (0  : y(6)  = n)  is  a closed  polyhedral  convex  set  for  any  tj  in 

T“(S).  Hence  the  ®(q)'s  form  a finite  cover  of  S.  Thus  if  Rule  2 is 

applied  for  some  t)  in  y(S)  infinitely  often,  then  0( q)  contains  a 
“ 1c  — 

cluster  point  0 of  the  © sequence.  Suppose  0 lies  in  the  3-face, 
then, in  fact,  Rule  2 would  have  terminated  the  algorithm  for  the  program 
(3, n)  the  first  time  q was  encountered,  since  0 yields  a quasi- 
equilibrium. □ 
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Obviously ^ if  K is  taken  as  the  set  of  integers  and  6 is 
large  or  if  v is  large^  implementations  of  Rule  2 is  prohibitively 
burdensome.  The  idea  is,  howeve^  to  determine  what  (3- face  contains 
in  its  interior  a cluster  point  of  the  0 and  then  to  run  the  test 
ever- so-of  ten  for  different  y(  0)  as  the  ©^  get  close  to  the  f3-face. 

Note  that  if  the  payments  vector  p(x?  nf  0)  for  each  0 in  0 
is  based  upon  the  same  x,  then  (x,  it)  is  an  equilibrium  where 
it  = £ jt^k.  Also,  since  the  sequence  Xto^-v  ©^  for  k = 2,  3,  ... 

contains  only  finitely  many  terms,  at  some  iteration  one  might  be  able 
to  guess  those  terms  which  are  in  X(©  ) infinitely  often.  In  this  case 
one  can  implement  a special  "tail"  routine  in  which  the  aggregate  linear 
program  and  its  dual  is  not  resolved. 


I 


8.  An  Example 

A preliminary  version  of  the  algorithm  of  Section  7 has  been  pro- 
grammed by  D.  Solow  and  used  to  compute  an  c-equilibrium  of  Mas-Colell': 
pure  trade  model  of  [9].  We  first  describe  the  model  and  then  the  com- 
putational results. 

The  trade  model  has  three  players  and  two  goods.  Each  trader 
has  one  unit  of  each  good  and  the  utility  functions  are,  respectively, 

min(x,2y) 

min(2x,y) 

min(4x,5y) 


where  x and  y are  the  quantities  of  the  two  goods  consumed. 

Note  that  equilibrium  prices  of  this  model  are  irrational.  Recasting 
the  problem  into  our  format  we  have  c^  = = c^  = (10000), 

bx  = b2  = b3  = (1  1 0 0 0 0 0 0) 


1 (0  1 0 0 o\ 

1 = 1 0 0 1 0 0 j 

.2  _ (l  -1  0 1 0\ 

A1  ^1  0-2  0 lj 


2 = fl  -2  0 1 0\ 

2 |l  0-1  0 lj 

.2  _ { 1-4  0 1 0^, 

3 U 0-5  0 1 I 


1 


and  each  A,  is  8x5. 


The  computational  results  are  summarized  in  the  following  table 


where  k is  the  iteration  number,  e is  that  of  the  e-equilibrium 


and  s-indicates  the  cumulative  number  of  simplex  pivots  required  to 


solve  the  aggregate  linear  program  (including  Phase  I for  initialization) 


Hence,  for  example,  a ( . 00008) -equilibrium  was  computed  with  56  simplex 
pivots  and  73  iterations. 

The  overall  run  time  was  about  four-tenths  of  a second  on  an 
IBM  370/168.  The  computation  behaved  essentially  as  expected  and  neither 
confirms  or  denies  that  the  algorithm  is  a practical  device  for  models 
in  the  target  area  of  five  traders  and  two  hundred  goods.  Applications 
to  more  serious  problems  will  be  reported  upon  at  a later  date. 

Observe  that  the  convergence  rate  is  linear  in  the  tail  wherein 
the  error  is  halfed  every  4 or  so  iterations.  If  the  vectors  which 
appear  in  u X(0  ) infinitely  often  can  be  identified,  then  quadratic 
convergence  can  be  obtained  by  using  Newton's  method  or,  for  example, 
Saigal  [23]. 
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